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ABSTRACT 

The directrix approach, which gives powerful insight to the 
geometric structure of solutions of the general Riccati equation, is 
developed. 

Burgers* Riccati equation is derived, and conclusions are drawn 
utilizing the directrix approach concerning the boundedness properties 
of solutions of this equation with certain restrictions on parameter 
values. 

Closed form solutions are developed for Burgers* Riccati equation 
for certain parameter values. A method is produced to obtain the verti¬ 
cal assymptote for unbounded solutions of Burgers* Riccati equation 
with certain restrictions on parameter values. Conclusions drawn from 
the application of the directrix method to Burgers* equation are veri¬ 
fied. 

A research bibliography for Riccati's Equation is Included. 
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I. INTRODUCTION TO BURGERS' RICCATl EQUATION 


A. PHYSICAL SETTING 

The Riccati non-linear differential equation has many applications 
in mathematical physics. In this thesis a relatively unknown application 
of the Riccati equation to turbulence is explored. 

Consider the equation 


V - t) 



( 1 ) 



where v =* v(y,t), U = U(t) and /JL and l) are parameters. Equation (1) 
with yM = 2 appeared in J. M. Burgers* theory of a model of turbulence 
(See (!]). In this model U(t) is the analogue of the mean motion of a 
fluid through a channel. The independent variable y represents the 
cross-channel coordinate and t represents time. If the domain of y is 
taken to be [o,b], then v(y,t) represents the secondary motion or tur¬ 
bulence of the fluid moving through the channel under some external 
force. 

If in Equation (1) U(t) = 0, then Equation (1) becomes 



( 2 ) 


Suppose that y€ [0,b], v(y,0) is some known function vo(y), and v(0,t) = 
v(b,t) = 0. In 1951 J. Cole (see [2]) published the general solution 
of this system in terras of the heat equation 


de - i) 

dt ^ 


where 


= -2-1) J_ 10 

G 
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With regard to this system, where y takes on values in a finite range, 
the theory of (2) becomes completely understood, since the theory of the 
heat equation is completely understood. 


If the eigenvalue-type restriction of the initial and final values 
of v(y,t) are lifted, and if y is allowed to take on values in an in¬ 
finite range, then v(y,t) represents the free turbulence in an infinite 
channel. Thus, the problem that Cole solved is actually a special case 
of the free turbulence problem. It is possible that an extension or 
generalization of Cole*s method of solution might lead to a complete 
solution of the free turbulence problem associated with equation (2). 

However, Equation (2) will be approached from another direction - 
via the geometry of the classical form of the Riccati equation. The 
approach is of such breadth as to apply to Riccati equations of greater 
generality than those associated with (2). 

To see how the connection with (2) and the Riccati equation is made, 
observe that in Ref. 1 a typical solution of (2) is put into the form 


Vf-l) 







where 



If the above transformation is applied to (2), a Riccati equation can be 
formed. 

B* MATHEMATICAL DERIVATION 

Thus, consider the following equation: 




(1) 
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with parameters yx , “i) . Let v( ) = V(‘^)/(t - to)^ where *] “ (y - yo)/ 

% 

(t - to) . Then, substituting above and using 

= - (nV' + V)/Z(i-t.)^ 

= VV (i~U) 

and V"/(i--tO^ 


where denotes ^ 


(1) becomes 

-n V'- y + uvv' = . 

a z ' 

Thus, 

f t)V'- V'^ + V 

^ 2 . z 

] =0 , 

and 

i)V'= i \/' - jV + 


2 

where C is 

an integration constant. 


In 

iiv' = i fv'-iV 

2 V /> 

• I'O 

let 

y - /A-l 
' Z-v) 


Then 

t) \J' = i) <j V M = 1 

- y v + y c' 


dY ^ ' 


or 

dV - _ 1 \/ + 

ciy 

2 c' . 

A" 

Now, since 

II 



dV - V' - yv + z c" 
dy 


Equation (4) will be referred to as Burgers' Riccati equation, or simply 
Burgers* equation. A more streamlined form in which Burgers' equation 


M 

t 





I 

I 


I 



will be handled is 


y'= , 

where C< and p are constants, and denotes ^ . 
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II. DEVELOPMENT OF DIRECTRIX APPROACH 


A. DEFINITION OF DIRECTRICES 

Consider the general Riccati equation 

where f(x), g(x), and h(x) are real-valued functions of x. For pur¬ 
poses of this thesis, h(x) is assumed to be non-zero for all x. Since 
the right side of (5) is quadratic in y, (5) can be factored into the 
equivalent form: 

r 

uCx')^ [-^(xH (gHx)-4f(x)^(x))']/2h(x) 
and i{y)~ [-g(x)-(g^(x)-4-p(x)j(x))'J/Ph(x). 

Let D be the set of all x such that u(x) and (x) are real-valued. For 

x^6D, and y(x) a solution of (6), the value of y*(x^) depends upon the 

position of y(x ) relative to u(x ) and (x ) (See Figure 1.). Clearly 
o o o 

y*(x ) =0 if and only if either y(x ) = u(x ) or y(x ) -i.(x ). Also, 
o o o o o 

y*(x ) ^ 0 if and only if J^(x ) y(x ) K. u(x ), and y*(x ) >“ 0 if 

o o o o o 

and only if either y(x^) > Knowing u(x) and 

^ (x), then, the solution may be "directed" at x^6D, in the sense that 
the slope is known at x^. Thus, u(x) and ^(x) "direct" the solution 
of (6) and will be referred to throughout this thesis as upper and lower 
directrices, respectively. 
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B. EFFECTS OF DIRECTRICES ON SOLUTIONS 


From an intuitive point of view u(x) appears to repel a solution 

y(x) in the sense that for x 6D, y(x ) > u(x ) implies that y'(x ) > 0, 

o oo o 

andi(x ) < y(x ) < u(x ) implies that y'(x ) ^ 0. Similarly, ^ (x) 
ooo o 

appears to attract a solution y(x) in the sense that y(x^)< implies 

that y'(x^) > 0» and, once again, (x^) < ^ implies that 

y'(x^) < 0. A simple and interesting case of the repelling property of 
u(x) and the attracting property of ^ (x) occurs when u(x) = C 2 >C^=^(x). 

Suppose now that u(x) = C^ > C^ a ^ (x), and that y(x) is a solution 
of (6). It is interesting to observe that y 2 (x) = C^ and yj^(x) a C^ are 
both solutions of (6), By a uniqueness argument, then, the only solution 
that can assume the value C^ is y 2 (x), and similarly for C^ and yj^(x). 

Thus, the solutions of (6) with u(x) a C^ > C^^h ^ (x) are partitioned 
into three distinct categories (not counting y^(x) and y 2 (x)). 

Suppose that y(x) takes on values greater than C^. Since y(x) ^ C^ 
implies that y'(x)> 0, as x-^ CX) , y(x) is repelled from u(x). In 
fact, since r , t \ *1 

f= ^)] , 


y(x) > C^ implies that y"(x) ^ 0, so that y(x) is concave upwards. 
Suppose that ^ (x) < y(x) < u(x). Then for y(x) > (C^^ + y' (x) ^ 0 

and y'‘(x) < 0. Thus, y(x) is repelled from u(x). For y (C^ + 

y*(x) < 0 and y"(x) > 0, so that y(x) is assymptotic to J (x) = C^ as. 
x->oo. That is, y(x) is attracted by (x). S uppose finally that 
y(x) < C^. Then y'(x) > 0 and y'‘(x) > 0 so that y(x) is assymptotic 
to J? (x) as x-^oo . Thus, in the case where u(x) s C^ > Cj^ s ^(x), 
solutions that are attracted to are bounded as x->oo and those 

solutions above u(x) are repelled and become infinite as x->oo, (See 
Figure 2). 
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If the directrices are allowed to vary with x, the bounded and un¬ 
bounded solutions of (6) can still be studied but are not as easily cate¬ 
gorized as in the constant*'directrix case. Specifically, if the directrices 
are such that solutions may intersect the directrices, perhaps even in¬ 
finitely many times, the regions in which solutions become unbounded are 
relatively more difficult to ascertain. Nevertheless, observations at 
specific values of x can always be made concerning the "direction" of the 
solution, and deep intuitive insight into the forms of solutions is gained 
using the geometric structure of the upper and lower directrices. 

C. APPLICATIONS 

1. Riccati "Lens " 

As an interesting and novel application of the directrix ap¬ 
proach, consider the problem of constructing directrices in order that 
the solution have certain properties. For example, it has already been 
shown that solutions can be repelled from or attracted to a certain value 
K^. In the first case, let u(x) = and (x) = where 
In the second case, let ^ (x) = and u(x) = where < K^. The 
degree of attraction or replusion can be adjusted by varying or K^. 

Thus, a sort of Riccati "lens" can be constructed, in that the solu¬ 
tions are either focused (attracted) or scattered (repelled). (See 
Figure 3). 

2. Oscillating Solutions 

Suppose now that an osciallating solution is desired. Con¬ 
sider (x) =* AsinBx + C and u(x) s K> A + C. Let y(x) be a solution 
of (6) and suppose that for some x^, j2 (x^) ^ Then 
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y*(Xo)<l 0* As X increases from x^, y(x) must eventually intersect 
Ji, (x) , since y*(x) 0 in this region. Since y(x) must intersect (x) 

at some point Xj^, where y*(xj^) = 0, J^'(Xj^) must be positive. As x in¬ 
creases from Xj^, y*(x) ^ 0, so that y(x) must necessarily intersect Jl (x) 
again, say at x^. But this puts y(x) for x > x^ between u(x) and J? (x) 
and the argument repeats. Furthermore, y(x) cannot assume a constant 
value for x > x since y*(x) = 0 implies that y(x) = K = u(x).for x > x. 
Thus, any solution y(x) < u(x) oscillates for x > x^. (See Figure 4). 

3. Random Directrices 

Another interesting problem is obtained by partitioning the 
x-axis into subintervals u^, for i » 1, 2, 3,..., and allowing u(x) and 
JJ, (x) to be defined for each subinterval. Let u(x) and (x) be assigned 
random constant values for each subinterval. (See Figure 5). The figure 
then shows but one realization of an ensemble of u and Jl segments. These 
U and J? values can have joint probability density functions of arbitrary 
type. Some interesting questions are: 

(a) Given an initial value for a solution y(x) of (2), 
what is the probability that y(x) will be unbounded as x-^oo ? 

(b) How does the answer to (a) depend upon the values that 
u(x) and j2.(x) can assume? 

4. Approximation of Continuous Directrices 

Consider, finally, two real functions fj^(x) and f 2 (x). Let 
the x-axis be partitioned as before, and define step function approxima¬ 
tions of fj^(x) and f 2 (x) on this partition. Let u(x) and XM be de¬ 
fined for each interval as the value of the step function approximations 
to f^(x) and f 2 (x), respectively. Note that given an initial value y^. 


14 










we can trace, through each subinterval, the solution of the Riccati 
equation defined by the values of u(x) and j2(x) on that subinterval. 
Suppose that the size of the subintervals approaches zero. Will this 
interval solution described above approach the solution through of the 
Riccati equation defined by and f 2 (x)? 

The consideration of all of these questions is beyond the scope of 
this thesis. They have been introduced to show the suggestive power of 
the directrix approach to the Riccati equation. 

5. ’ ’Inverse" Problem 

It is widely known of the Riccati equation that a general solu¬ 
tion can be obtained from a non-trivial particular solution by means of 
two quadratures. Thus, the problem of finding general solutions reduces 
to that of finding non-trivial particular solutions. In some cases a 
simple regrouping of the terras in a Riccati equation can lead to a parti¬ 
cular solution by inspection. Sugai [Ref. 3], for example, gives some 
particular solutions for cases in which the coefficients in Riccati*s 
equation satisfy a prescribed interrelationship. The factoring of 
Riccati*s equation, however, lends further insight to the problem of 
finding particular solutions. Suppose that f(x) is a real-valued function. 
A Riccati equation that is satisfied by f(x) can be found by making cer¬ 
tain assumptions about the directrices of the Riccati equation. Specifi¬ 
cally, if one directrix is required to be a certain function, then the 
other directrix can be found through a computation involving f(x) and the 
known directrix. Then, since both directrices are specified, the Riccati 
equation satisfied by f(x) is completely specified. 

The directrix approach, then, suggests the following "inverse" prob¬ 
lem. Given one directrix, d^(x), and a function f(x), what must d 2 (x) 
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be in order that f(x) satisfy the Riccati equation y* =» (y-d^) (y-d^) ? 

Here the notation u(x) and ^(x) for directrices has been changed to 
dj^(x) and ^2^^^ since the "upper" and "lower" ordering may be lost during 
crossings, i.e., u(x) (or j2.(5c)) may turn out to be defined as d^(x) in 
some region and as d 2 (x) in another region. Solving for d^, 

da f''- + dif , 

C1.--P 

Consider, for example, the Riccati equation 

' = - da) (8) 

That is, d^ ~ 0 and, from (7), 

da = f - -py-f 

Suppose, for example, that f(x) is a polynomial of degree n > 0 such 

that f(x) ^ 0 when x€D' for sane domain D*. What must d^Cx) be in 

order that f(x) satisfy a Riccati equation of type (8)? From (8) it is 

clear that deg(f-d^) = Since ^ must be an integer, n =* 1. 

z n n 

Hence, f d^ * C. But f - d^ =* C implies that f* * Cf; that is, 
f * exp(Cx). Thus (8) has no polynomial solutions of degree n > 0. 

This example serves actually two purposes. The first purpose, of 
course, is to point out that, using this method to pair particular solu¬ 
tions to equations, not every function can be realized as a particular 
solution of a Riccati equation of specified form. The second purpose is 
that the example leads into the problem of finding polynomial particular 
solutions of a Riccati equation whose coefficients are polynomials. 

This problem was solved :n 1954 by Campbell and Coulomb [Ref. 4], for 
the Riccati equation 
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( 9 ) 


A = Bo + Bi ^ 

where A, B^ are polynomials in x. In Ref* 4 criteria are pre¬ 

sented for the existence of, and an algorithm is presented for finding 
all polynomial solutions of (9). For the solution to the problem for 
A H B=l, see Rainville [Ref. 5). 

6. Tables of Particular Solutions 

Suppose, now, that f(x) == a/x and d^ = 0. From (7), d^ * 
(l+a)/x so that f(x) satisfies 

Using this procedure, with further choices of f(x) and d^(x), a table 
of specific Riccati equations and the corresponding particular solutions 
could be constructed. Thus, the directrix approach provides a systematic 
procedure for constructing the Riccati equations associated with parti¬ 
cular choices of one of the directrices of the equation and an arbitrary 
solution. 

D. EXAMPLE 

1. Template Solution 

Suppose, on the other hand, that a specific Riccati equation 
is given, and that the solutions are to be constructed. The directrix 
approach yields insight to the behavior of solutions in certain regions 
of the plane. It is instructive to consider the example that follows. 
Consider the Riccati equation 

^ + b ^ + C 

where b and c are real constants. 



17 






i 




An equivalent fonn of (10) is 

^'= 

where ^ 

r, = [ - b + ( b - 4c)‘ ] / 2 , 

r. = [-b-(b‘-4c)*]/2 . 

Suppose that b - 4c 0. Then, since the directrices are identically 

constant, the solutions of (1) can be categorized as in Figure 2. It is 

interesting to observe that y'(x) does not depend upon x. Thus, for any 

point (x ,y ), the unique solution that passes through the point (x ,y ) 
o o o o 

must have the same slope as the unique solution that passes through the 

point (x,y^) for any x. That is, given a solution of (10), any other 

solution of (10) in the same region can be realized as for some 

X . A "standard" curve can thus be constructed for each of the three 
o 

regions, whose shape is determined only by the constants b and c in (6). 
This "standard" curve or "template" can be translated along the x-axis 
to represent the unique solution of (10) passing through any particular 
point in the region of definition of the "standard" curve. (See Figure 


6 .) 

2. Verification of "Template" Solutions 

The closed form solution of (10) can be found, by direct inte¬ 
gration of (10), to be 


where 


u(x) = ' T, - ft C u (y) 

^ l-Cu(x) 

u.(x) = exp[ (r.-Tz) X ] , 

X = X- Xo , 
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and 


Notice that when Cu(x) 
when 


c = i^(x.)- r, 

1, y(x) is infinite. That is, y(x) is finite 


exp [ (r,- rj X ] = _1 


C 

For a solution y(x) > 0 4 C 4 1, so that 1 ^ 1/C 4 CXD 

Since “ ^2 ^ there is an x such that exp[ (r^^-r^)?^] = ^ , so that 

y(x) is assymptotic to x = x as x xr Similarly, for a solution y(x) 

^ r^, 1 4. C 4 00 ^ so that 0 4 ~ 4 1, In this case x ^ 0 implies 

that exp((r,-r )x) y 1 ^ Thus, y(x) is not infinite for any x ^ 0. 
L Z O 

A — A _ 

However, there is an x 4 0, such that x = x implies that exp((rj^-r 2 )x) = 

1 + — A 

That is, as X X , y(x) becomes asymptotic to x == x. 

c 

Furthermore, for a solution y(x) such that 4 y(x) 4 ^ 

simple transformation of variables reduces (11) into a familiar form. 
Since any one solution in this region is representative of all of them, 
let y(x^) = (r^+r^)/2. Let y(x) * y(x) - (rj^+r2)/2 where x = x - x^. 


Then 


u(x)= - a tanh(ax) 

0 2 p 


( 12 ) 


where a = r^^ - r^. 


The "standard" curves or "templates" can be expressed for the three 
distinct regions as special forms of Equation (11). For example, letting 
C * 1, the "standard" curves obtained are, 

u(x) = fi - n exp(ax ) 

wherea-r^.r^ ^ l-eXp(dX') 

for U(x) ^ [ n ) Til ) 
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and 



where a * 

for n. ^ i^cx) L r, 

These "standard” solutions are exactly those represented in Fig. 6. 

2 

Consider now the possibility that b - 4c 4 0. That is, suppose 
that y* = (y-rj^) (y-r^) with r^^ = c + di and r^ = c - di. Direct inte¬ 
gration yields 


<^(x)= d tan[d(x+K)] + c 


(13) 


where K is an integration constant. Consider equation (12) with a = r^^ 
r^ ** 2di. Then (12) becomes 

i^(x)- ^ - -di tanh[di(x-Xo)] 

= ’-d tan[-d(x-Xo)] 

Finally, l^(x)= c\ t dH (dX “ clXo) % . (13*) 

Note that (13*) agrees with (13) with-K - x and c = Equation (13) 

o / 

is graphed in Fig. 7 for K =* 0. Note that (13) resembles a 90® counter 

clockwise rotation of (12). 

2 

For the case where b - 4c ^ 0, then, the directrix approach gave 
insight into certain properties of solutions of (10). Specifically, 
boundedness could be discussed in terms of the region in which the solu 
tion appeared. These regions were determined by the directrices. The 
independence of y' with respec t to x sugf.asted a "template" representa¬ 
tion of the solutions of (10). The directrix approach, therefore, can 
give strong geometrical insight to the solutions of (10). What insight 
can this approach give to Burgers' Riccati equation? 
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III. APPLICATION OF DIRECTRIX APPROACH TO 

burgkrs* riccati equation 


A. CASE I: y3 = 0 

Now, consider the Burgers *-Riccati equation 


from the directrix point of view. The case where ** 0 will be 
considered first. 

Consider the equation 

^ - 2.o(x . (15) 

An equivalent form of (15) is 


^ ^ 


+ yS ') o(> 0 (1^) 




doc X 



i^( Z(xx) 


(16) 


Note that y * 0 satisfies (16) so that no non-zero solution of (16) can 
cross the x-axis. Hence, only solutions y(x) such that y(x) > 0 for 
all X are considered in the following arguments. 

The directrices for (16) are 

/ N i 2.(xx 

u(x')= I Q 

and 


defined by 

for X > O 
for O 


i(x) = 



for X O 
for x<0 


For y(x) > 0 a solution of (16), y*(x) is described as follows: 
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a 


'(x) 


= 0 
< 0 
> 0 


for 


for 

0< 1 

for 

^(x) 


or u- 2ax 


> 2 c( X 


For r a real number, define 

= (^x,Lj)ly>0 and Eo(x^-r = o]- . 

Then, the solution y(x) of (16) that passes through (x,y) 6 must 

have slope equal to r when x = x. Notice that the equation 

^^-2cxxi^“r= 0 

can be written as 

[. ij- ojCx)] [<^-v(x)] = 0 

uj(xl = cxx + (o(^x^ + rV 

I 

v(x)= OCX - (o('x^ + r)" 


where 


and 


(17) 


In particular,= 2oCx. The family of curves defined by the relation 
as r ranges over the real numbers is represented in Figure 8. 

Notice that for r > 0, defines a single-valued function for posi¬ 
tive values of y. For r > 0 let ^^(x) denote that function, -^^(x) 
is in fact equal to u(x) for positive y and r > 0. The slope of -^^(x) 
is given by: 


4r(x)= cx + X 

(x^x^ + r)'"" 


(18) 
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Note that 04 ^ 2<X for all x, and that rx 0. 

Also, -^^(x), given by 


Note that 0 4 


is positive for all x. 




Suppose y(x) is a solution of (16) such that for some x > 0, y(x ) 

o o 

4 2 (X x^. Then ^ 0, and as x-^ oo , y(x)—> 0; that is, y(x) 

is bounded as x—> oo . Thus, any solution y(x) such that y(x ) 4. 
for some x^ > 0 is bounded as x —^ cc> , Moreover, since no bounded solu¬ 
tion can be entirely above y = 2 OC x for all x > 0, then every bounded 

solution must lie between the directrices for all x x for some x . 

o o 

That is, the existence of an x^ such that x > x^ implies that y(x) —=> 0 
as X— >oo completely characterizes bounded solutions of (16). 


Now, suppose that for some x^ > 0, 2 DC x^ 4 y(x^) ^ 

Then 0 4 y*(x^) ^ OC , From equation (16), y"(x) is found to be 


t^"(x) = 2 [ t^( 1 ^'- (x) - 


(19) 


Thus, y*’(x^) 4 0 since y* - OC. — 0 and - OCxy* 4 0. Hence, y(x) must 
intersect 2 oC x for some x > 0. By a previous argument y(x) must then 
be bounded as x —> co , 


Suppose that for some x > 0, y(x)^ y*(x) ^ 2a. 


Since cX 4. 2a, y‘(x) must be increasing as x increases from 

X. Alternately, from (19), 



Thus, y(x) cannot intersect y = 2 x and is unbounded as x—> co 
To summarize the results (See Figure 9): 
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(1) If for sorae > 0, y(x^) ~ then y(x) is 

bounded as x —> oo . 

(2) If for some x > 0, y(x) ^ 
bounded as x —>> co . 

B. CASE II; ^ =• 2C< 

It is interesting to notice that the solutions of Burgers' Riccati 
equation are symmetric with respect to the origin. That is, 

4iLl) =: (-yf- £ci(.(-x)(-M)-(-4 = u^-2(X.xu-ts = 4^ 
cl(-x) ^ ^ ^ ^ ^ dX 

It is thus sufficient to consider only those solutions in the half-plane 
where x ^ 0. There is, by symmetry, a solution in this half-plane which 
corresponds to any solution entirely contained in the half-plane deter¬ 
mined by X <. 0. This fact is helpful in considering Equation (14) with 
/Sjt 0. 

Equation (14) will now be considered for ^ =* 20^ , that is 

+ E<x . (20) 

As in the case where * 0, define a family of curves; 

I 1 ^otxy + 2(X-r = 0^ 

Notice that for y ^ 0, (B^ = ^2 • fact, for y > 0, the family 

of curves is merely the fa tiily of curves with the subscripts in¬ 
creased by 2oc . In particular, the directrices of (20) are actually 
Figure 9.) 

Notice that y = 0 does not satisfy (20). Thus, solutions may cross 
the x-axis, and in fact, cross with slope equal to 2 . However, y « Hc(X 

does satisfy (20) so that no solution can cross the line y ® 2oCx, Notice 
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also that y =» 2 (X x is precisely (x), since the solution y =» 2a x 

2 a 

has slope equal to 2^ for all x. The following discussion pertains to 
all solutions y(x) of (20) for which y(x) > 0 for all x. 

Suppose that for some x ^0, y(x) = 0. Then for some x' > x, y*(x') 
0; that is, y(x') 6 (B^. If not, then y^ 2cxx for x > x, which contra¬ 
dicts that y(x) = 0. Thus, since x > x* forces y'(x) ^ 0, y(x) must be 
bounded as x->>co . Notice also that if y(x) = 0 for x ^ 0, then y*(x)> 
for X > X. Since for x > x, y(x) > 2a x, then y(x) must be un¬ 
bounded as X— >oo. 

The boundedness of a solution y(x) of (20) as x —^ oo can thus be 
determined by the sign of the x-intercept of y(x). Notice that every 
solution y(x) has an x-intercept so that this characterization for 
boundedness of solutions is meaningful. Moreover, note that if y(x) 
is bounded as x —> co , it is unbounded as x —> ^^9 sod vice-versa 
(that is, interchanging 00 and -oD), except for the case y = 2ax. 
Furthermore, the characterization for bounded solutions could also be 
given in terms of the y-intercept. That is, y(0) > 0 implies that y(x) 
must have a negative x-intercept, so that y(0) 0 leads to an un¬ 

bounded solution as x-^oo. Similarly, y(0) ^ 0 leads to a bounded 
solution as x—><». 

To summarize the results (See Figure 11): 

(1) If y(0) ^ 0, y(x) is bounded as x-^oo, 

(2) If y(0) > 0, y(x) is unbounded as x->oo. 

C. CASE III: yS>0 

If were allowed to take on all values between 2c< and zero, that 
point p on the y-axis corresponding to =* 2c< , initial values above 
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which cause the solution to be unbounded and below which cause the solu¬ 


tion to be bounded as x—>oo, should take on a continuum of values be¬ 


tween zero and the corresponding point q for ^ = 0. It is known that 
q lies in the interval ( , \f2u ), but its exact value is not 


known through this geometric approach. As jS takes on values above 2CX , 
p moves down the y-axis until eventually p =* - OO . 

Note that increasing /3 causes a shift in the labeling of the family 



for r €: IR . In particular, the lines y = 2 x and 



y = 0 are always 


It should be expected, then, for 


some value of ^ 2 , that every solution which passes through the 

y-axis will be unbounded. The particular symmetry properties of the 
solutions, at least in the present case, suggests that if a solution is 
bounded for x-^oo and /S>2oc , then it is unbounded as x- oo , This 

statement holds, of course, for , since if y(x) is bounded as 

x-»oo , then y(0) 4 0. Thus -y(0) > 0, which implies that -y(x) is 

unbounded as x—^ - OO . For ^ 2oC , p 4 0, so that for any bounded 

solution y(x) as x~>oo , y(0) 4 p. Thus, -y(0)>-p, so that -y(x) must 
be unbounded as x —^ - oo . Note that the converse of the result is 

not true in general. That is, some solutions are unbounded both as x-^co 
and also as x —^ oo . A simple example of this for = 2^ is the 
particular solution y =* 2 OC x. 
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IV. CLOSED FORM SOLUTIONS 


A. CASE I: 0 

The closed form solution to Burgers' Riccati equation can be obtained 
with certain restrictions on o( and yS . 

Suppose = 0. That is, a solution of 


„5, 

is desired. The closed form solution is known (see Kamke [Ref. 6] and can 
be obtained from 

-L ex-pL io 2oCtdt] io expf-io ScCScisldt 

fx) 

Thus, 

_ exp [- So 2ixtcltl 

r ’ ~ —TTt-1 

Jo exp[- So 2<xsds] cl t 


B, CASE II: /S=: 

Suppose that/S ® 2 c< . It was observed that y = 2 c< x is a parti¬ 
cular solution of 


1 / — 


} = ^ 


- 2(xxi^ 4- 2(X 


( 20 ) 


Let y ^ 2o<,x. Then, substituting into (16), 

a' - 2 ccx a - 1 = o 


The solution of this non-homogeneous linear first order equation is 


given by 


u(x) = exp[- So dotldt"] [u(o) - $0 exp [ L dxsclsl d tj 
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Since 


JL 

U. 


■t Z.O(,K 


) 





exT> 


(o 2t<tdt] 

1 _ 1 - 2cXX . 

u(o)- S 

0 exp[ 

loZo^ScIsldt (22) 


Equation (22) is the closed form solution of (20). 

For other values of ^ , the solutions of Burgers' Riccati equation 
are not as easily found. For certain classes of yS , however, solutions 
can be found in terras of known functions. These solutions will be 
developed in the following pages. 

C. METHOD OF FROBENIUS 

It is well known that the Riccati equation is equivalent to the 
linear second order differential equation. That is, to each Riccati 
equation there corresponds a linear differential equation of second 
order, and vice-versa. In seeking solutions to 

= - E(xx 

it is advantageous to study the second order linear form of (14). If, 
in particular, y “ -u*/u, then (14) becomes 

u" + ZoCX u' + /3 U •= 0 , (23) 

The method of Frobenius (e.g., Wylie [Ref. 7]0 can be applied to 
(23). For example, suppose that 

u = X*" [ do + a.X + + • • • ] (24) 




(14) 
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Substituting for u back in (23) and equating coefficients of like powers 
of X, the coefficients a^, i = 1,2,3,..., can be found. The roots of the 
indicial equation turn out to be r = 0 and r = 1. Since the two roots 
differ by an integer, only one series solution can result. Thus, r - 0 
is chosen for simplicity. The values a^ and a^ are arbitrary constants, 
and the values of the other coefficients are 

= ( -1) do [2K(Z.K)-f/S][?o((2K-^)+/3]--- 

( + 1 )! 


for K = 0,1,2,3,... and 

(-1)' di [^P(( 2 K-i)+p] K- 3) H-^]- [ gg 

(2K+ 1)! 

for K * 1,2,3,... 

Thus, one solution of (23) is (24) with the ^^^*8 defined as above. To 
obtain a second solution, assume that u(x) * (|) (x) Uj^(x), where Uj^(x) 
is the solution of (23) corresponding to (24) with r = 0. Substituting 
for u(x) in (23), and noting that u^(x) satisfies (23), (|) (x) must then 

satisfy / \ j , 

u, $ + ( 2 Ui' + ^c< X Ui) (|) - 0 . 

<i>w = c5: . K. 

Then, u *(t)uj^, so that the solution of (14) y ** -(u‘)/u becomes. 




ul(x) __ C exp (-cxx^) 


• (25) 
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Thus, (25) is the general solution of (14). 


It is interesting to note that for /5 » 0, all of the even indexed 
coefficients of (24) become zero. If a^^ is taken to be zero, u^^ s a^. 
Thus, (25) reduces to 

u(x) = - ex? (-Kx^) _ 

^ (K/doC) H- ^0 ex'p(-0(i‘^)clt 

which agrees with (21) if the integration constant K is taken to be 
zero. 

Similarly, for y6 = 2iX , if a^^ is taken to be zero, = exp(-o(x^). 
Substituting for u^ in (25) gives 

u(x) zi Ecxx + e><T>M) 

^ (-k/C)- exp(o(,i^)dt 

which agrees with (22) where u(0) = — . 


D o ABDELKADER'S METHOD 

The closed form solution can be obtained for an infinite class of 
*8. The following method is due to Abdelkader [Ref. 8], Consider 

'i' = • (u> 

Let y(x) = a(x) - V*(x)/ \r(x). Substituting into (10) for y, 

v" - [ 2a - 2«.xlu' + [^ + o<.(-E«x)+ 4c<^x^ - a'lv = 0, 

or. 


v'' - A(x)v' + B(x)ir = 0 


( 26 ) 


If a(x) can be chosen so that A'(x) - B(x), then (26) can be solved 
easily for 'VJ(x). That is, (26) becomes 

v" - [ A(x) vl ^ = 0 
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from which 


\r' = A(x)ir -i- Cl 

can be solved easily for \r(x). Then the general solution of (14) can 
be written 

(x) = 9(X) - A(x) - ^ 




v(x) 


The function a(x) is called the "conjugate" to y(x) and satisfies 

d' = ~ di^ -h ScKXd - ( /3- Eoc) , (27) 

If /S = 2<x , (27) can be solved easily for a(x) (see Kamke [Ref. 6]). 

If ^ = 0, then a(x) = 2 oc x is a particular solution of (27), which is 
all that is needed. 

If (27) cannot be solved easily as is the case if = 0 or = 2(X , 
then to repeat the above procedure on equation (27) would merely give 
y(x) as the "conjugate" to a(x). The transform-?.tion UJ = 1/a yields a 
distinct Riccati equation, namely 


UJ 


^ = 1 - 2o(,XUr + (/3 -£c<)lo‘ 


(28) 


Let (a) = a^^ - u'/u so that (28) goes into 


U 


U 


- C~2c)CX+ E 9, (/5--2e<.)3LL' -I- (jS- 1 + + ]u~0 


- A,(x)u' + Bi(x)u = o 


If aj(x) can be chosen so that A*(x) = B(x), then u(x) can be found as 
v(x) was found in (26). Then a^(x), the "conjugate" of (x), must 


satisfy 


a.' = -(/5-2oi)af + Ec<xa. - . (25) 

(/3 ~ 2 oc) 
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Notice that it ^ ^ 4cX. , (29) can be easily solved and U)(x) can be 
written 

U)(x) = d,(x)- Ai(x) - ^ 

u(x) 

Tracing back, y(x) can be found. Note that can be taken to be zero, 
since only particular solutions are needed at all except the last step. 

It is interesting to note the form of succeeding conjugates a^(x). 
For example, a^Cx) turns out to satisfy the equation 

a/= _(All?0 3.1 + ^0() 

(/3-4o() 

In general, if a^(x) satisfies 

a( = - Aaf + ^axa; - B , 

then a^^^(x) satisfies 

a;l, = -B a;', + 2otxa«, - A 

B 

where C - - 2(i+l)C<. ). Thus, it ^ - 2ncX , then the (n+l)st con¬ 

jugate, a^(x), can be solved easily. Working backwards, y(x) is event¬ 
ually obtained. 

Thus, this infinite class of jQ 's leads to a closed form solu 

tion of (14) and is defined by ^ = 2noC for n a natural number. 

E. HERMITE POLYNOMIAL SOLUTIONS 
Consider, now, the equation 

= (14) 

and apply the transformation y = / u;(t). t = \fp X. 

Then, (10) becomes 

- UJ^+2.^too+l , ojKere X = fL . 
dt /3 
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Thus, the number of parameters has been reduced from two to one. Let 
10 = -V'/V to obtain the second order form 


V"- 2 )St V' + V = O , "here ' = ^ 

Now, letting "t = % 

^V"- E W' + V = 0 

where differentiation now is with respect to 'S . Dividing by ^ , 

V" - 2 i V + 1 V = O . (3 

Now observe that the equation 




0 j for n a natural number, 

has H (x), the Hermite polynomial order n, as a particular solution, 
n 

Thus, for ( 2 ) ^ = 2n for n a natural number, equation (30) has\/( t)= 
Hn( ^ ) as a particular solution. The condition ( 2 ) ^ * 2n can be re¬ 
written as /^ / OC * 2n or =* 2n(\ . Thus, when yfi = 2n0< , for n a 
natural number, equation (14) has a solution in terms of Hermite poly¬ 
nomials. Notice that the condition that ^8 - 2n«, is precisely the 
restriction upon the class of solvable Burgers' equations using 
Abdelkader's method. 


Fc HYPERGEOMETRIC FUNCTION SOLUTIONS 

Still another approach to Burgers' equation comes from its linear 
second order form. Letting y = -u'/u, 

U." 4- ^CXXU' + yS U = 0 , (31) 
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Equation (31) has the solution (see Kamke, [Ref. 6], and [Ref, 9]). 

u(x) = X ‘ exp ) M (1- ^ , c<x^ ) 

z 

where M(a,b,x), the Kummer's Function, satisfies 

Xw" + ( b-x)tAj' - auj = o 

Burgers' equation thus has solutions in terms of confluent bypergeo- 
metric functions. 
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V. CONCLUSIONS 


A, BOUNDEDNESS PROPERTIES FOR = 2n 

It is interesting to note that if, in Burgers’ equation, y3 = 2n(X 

then a particular solution, y., is known in terras of H (x). Letting 

1 n 

y * + 1/v in (10), v raust satisfy 


\r' + (-Zc^x+2.i^j)u = -1 

\;(x) = C — ^0 A(’t)clt 
A(x) 

exp [ L (~2.oct + £y|(i)) dt] . 

Then the solution of (14) becoraes 


That is, 


whe re 




1 + . 


A(x) 


CA(t)clt 


(32) 


From (32), y(0) = y^(0) +~ and 


0= ( i^(o)- 


Note that y^(0) is known. Also note from (32) that y(x) becomes un¬ 
bounded when 


c= A(t)dt 

That is, given y(0) the constant C is completely specified. Then the 
solution y(x) corresponding to y(0) will become infinite for that value 
of X for which 


1 


Jo A(-t)clt 


(33) 









Therefore, the evaluation of the integral 


l!A(t)dt 


will yield the 


value of X at which an unbounded solution has a vertical assymptote. If 
no such value x exists as x —^ <x> , then y(x) must be a bounded solution 
Note that formula (33) is valid for any particular solution of (14) 
In particular, for ^ = 2n^ , formula (33) will tell if y(x) is un¬ 
bounded for a particular y(0) and in fact will give the assymptotic 
value for y(x) in the case that y(x) is unbounded. 

In this way the geometric directrix approach to the Riccati equa¬ 
tion is supplemented. That approach, it will be recalled, is applicable 
whenever a moving point has a pair of directrix points above, below or 
around it to guide it. However, the directrix will not be present if 
the roots of the factored Riccati are complex. In this case the solu¬ 
tion becomes unbounded. It is precisely in this case, that the above 
information will be useful. 

B. VERIFICATION OF GEOMETRICAL OBSERVATIONS 

The geometrical observations made for the cases ^ = 0 and ^6 = 2C< 
will now be evaluated in light of formula (33). 


For yS =* 0, it was observed that y(0) 4. 



implied that y(x) 


was bounded as x—> oo . On the other 



meant that 


y(x) was unbounded as x 



s X —> CO . Thus, for the /3 = 0 case, y^(x) = 0 
-2 c< t dt). From (33) j 


Notice that for s = 



1 
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Thus, 


fX 1 vpn^ 

Q I i A(Odt ^ rH 

Hence, y(x) will be unbounded and have a vertical assymptote if and only 

if 2 ^TT/oC ^ y(0) 4 CO ; that is, if and only if vF^E* ^ y(0) 4 oO . 

_ _ V TT 

Notice that 4 \f, so that y(0) ^ implies that y(x) 

^ TT 

has no vertical assymptote as x —^ oo , and, hence, must be bounded. 
Moreover, , so that y(0) > implies that y(x) 

has a vertical assymptote. Thus, the analytic expression has produced a 
sharp bound for the initial values of solutions that are bounded as x-^ co 
The geometric insight gained from the directrix approach, however, can be 
intuitively helpful not only to determine boundedness but also to visualize 
the general behavior of solutions as the initial value varies. The con¬ 
clusions drawn from the geometric approach, though not as precise as the 
analytic in this case, are nevertheless verified. 

For ^ = 2c< it was observed that any positive y(0) led to an un¬ 
bounded solution, and that any negative y(0) led to a bounded solution. 
Here, unlike the case where ^ 0, a sharp bound was obtained from the 

directrix approach, namely the value y(0) * 0. In this case, y^ =* 2 OC x 
and A(x) = exp ( 4-= exp( x^) . Then, from (33), since 

-L = Srexf(«-tOdt 

S'”' , 

2 

Notice that 0 4 exp( 0(. t )dt 4 . Thus, any positive y(0) will 

have a vertical assymptote and any negative y(0) will not have vertical 
assymptote, as x —> co . Hence, the observations are again valid. 

Notice that in this case the directrix approach led to the same sharp 
bound as was produced in the analytic argument. 


37 



Figure 1 - Position of y(x ) relative to l(x ) and u(x ) 

o ' o o 


y 
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Figure 2 - Constant directrices u(x) = C^; l(x) = 


y 
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Figure 3 - Riccati ’’lens 


A =» Focusing Region 
B =» Scattering Region 


y 
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Figure 4 - Oscillating Solution 
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Figure 5 - Random constant directrices 
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Figure 6 - "Template" solution through 


y 



43 











Figure 7 - Complex directrices u(x) = C + di; l(x) = C-di 



44 



















































Figure 8 - Family of 


curves 
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Figure 9 - Examples of bounded and unbounded solutions for 

/S=o 
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Figure 10 - Family of curves 
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Figure 11 - Examples of bounded and unbounded solutions 


as X—> CO for yS « 2 (X 
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